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Abstract 

We present an analysis of charm quark fragmentation at 10.6 GeV, based on a data sample of 
103 fb _1 collected by the Belle detector at the KEKB accelerator. We consider fragmentation 
into the main charmed hadron ground states, namely D°, D + , Df and A+, as well as the 
excited states D*° and D* + . The fragmentation functions are important to measure as they 
describe processes at a low energy scale, where calculations in perturbation theory lead to large 
uncertainties. Fragmentation functions can also be used as input distributions for Monte Carlo 
generators. Additionally, we determine the average number of these charmed hadrons produced 
per B decay at the T(4S) resonance and measure the distribution of their production angle in 
e + e~ annihilation events and in B decays. 

PACS numbers: 13.66.Bc, 13.66.Jn, 14.20.Lq, 14.40.Lb 



4 



I. INTRODUCTION 



Over recent years perturbative quantum chromodynamics (pQCD) has shown impres- 
sive agreement with various inclusive measurements at e + e~ colliders at many center-of- 
mass energies (CME) ranging from 14 GeV up to 206 GeV. These measurements utilised 
variables called event shapes or jet rates, see |l| for such an analysis. These are inclusive 
variables, whose values are calculated from the four-momenta of all particles in an event. 

Other properties, such as the momentum spectra of charged or neutral particles, have 
also been measured, but their prediction has proven to be more difficult. The necessary 
calculations have to cover the entire energy range from the production of the partons 
at the CME down to the scale of the hadron masses (typically 1 GeV/c 2 ), at which 
hadronisation occurs. Typically, powers of the form log (Q 2 /m 2 ) arise when quark masses 
are taken into account, making pQCD calculations difficult to interpret. 

Attempts have been made to extend the applicable range of pQCD to lower scales. 
These attempts have to be validated, for example by comparing so-called fragmentation 
functions. Due to the scaling violation of QCD, a fundamental property of this theory, 
the fragmentation function for a given particle depends explicitly on the CME. This 
energy dependence must follow the Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) 
(2I evolution equations. 

Thus, the fragmentation functions have to be properly evolved. Monte Carlo (MC) 
generators which include this scaling can be used instead of analytical evolution. Common 
MC generators which include this scaling are JETSET (its variant) PYTHIA j^J and 
HERWIG 0. 

These MC generators are also needed to model hadronisation, the transition of par- 
tons into hadrons, which cannot be calculated from first principles within QCD. Various 
models are implemented in MC generators. These can be distinguished by comparing the 
(identified) heavy hadron momentum spectra predicted by each model to the spectra seen 



5 



in data. 

Fragmentation functions for heavy quarks are attractive both experimentally and the- 
oretically. Concerning theory, mass effects in the matrix elements only have to be con- 
sidered for the heavy quark; in the limit of mu g ht — > 0, a pQCD calculation based on an 
effective Lagrangian reduces the complexity of the calculation compared to the case of 
light quark fragmentation. 

Experimentally, it is important to measure heavy quark fragmentation functions as 
their shapes are different from the corresponding functions for light quarks; such a mea- 
surement is furthermore straightforward, as very often hadrons containing heavy quarks 
can easily be identified. Since the production of heavy quarks is strongly suppressed 
in both the perturbative splitting of one parton to many partons (the so-called "par- 
ton shower") and in hadronisation, a heavy quark found in an event will most likely be 
produced in the primary interaction. 

At LEP and SLD, b quark fragmentation functions have been measured with high preci- 
sion f| Q, 0, E| . These measurements found that these fragmentation functions are in fact 
close to the ones of light quarks, suggesting that one combined model for all five flavours 
might describe the measured momentum spectra better than functions which have been 
introduced for heavy quarks alone. These collaborations have also published measure- 



ments of c quark fragmentation functions 



bj, but with large statistical uncvrtaintirs 



due to the small product of the branching fraction and reconstruction efficiency for the 
various final states. Some commonly used fragmentation functions are described by the 
models of Peterson et al. , of Kartvelishvili et al. [3] and of Collins and Spillerjb^ , 



as well as by the models of the Lund group j3| and one of its variants by Bowler 



1( 



For charm quark fragmentation functions at lower energies, the most recent published 
results for D°, D + , D*° and D* + are those of CLEO The analysis presented here has 
better statistical precision as the data sample is five times larger. Other measurements 
are more than 10 years old ; their data sample is over three orders of magnitude 
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smaller than that used in this analysis. The systematic uncertainties are reduced signifi- 
cantly and are comparable to those in ^?|. For a recent review of fragmentation function 
measurements and theory, see 20]. 

A measurement of D° and D + performed by the same experiment on the same data 
set allows for an easy comparison of charged meson production rates and momentum 
spectra, as well as a comparison of the momentum- dependent production of secondary- 
to primary-produced mesons. The measurement of the excited states D*° and D* + allows 
the determination of the feed-down contribution to the ground states D° and D + and also 
a momentum- dependent determination of V/(V + P), the ratio of the production rates 
of vector and the sum of vector and pseudo-scalar mesons. A comparison between 
production, and the production of D° and D + , can be used to determine the fraction of 
s quark production in hadronisation. Comparing the results for the A+ to those of the D 
mesons makes a study of the baryon production mechanism possible. 

In addition to charm fragmentation in the e + e~— >cc continuum, charmed hadrons 
in e + e~ annihilation events can be produced in decays of b-hadrons. The dataset for 
this analysis includes events above the production threshold for BB pairs, at the T(4S) 
resonance, so the lower momentum hadrons include contributions from B° and B + decays. 
This allows a measurement of the production rate of charmed hadrons in B-meson decay. 

II. DATA SAMPLE AND EVENT SELECTION 

This analysis uses data recorded at the Belle detector at the KEKB accelerator. The 
KEKB e + e _ collider is a pair of storage rings for electrons and positrons with asymmetric 
energies, 8.0 GeV (e~) and 3.5 GeV (e + ), and a single intersection point with a 22 mrad 
crossing angle. The beam energies are tuned to produce an available CME of \fs = 
10.58 GeV, corresponding to the mass of the T(4S). A detailed description can be found 



in 



2l|. 
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The Belle detector covers a solid angle of almost Att. Closest to the interaction point 
is a high resolution silicon micro- vertex detector (SVD). It is surrounded by the central 
drift chamber (CDC). Two dedicated particle identification systems, the aerogel Cerenkov 
counter (ACC) and the time-of-flight system (TOF), are mounted between the CDC and 
the CsI(Tl) crystal electromagnetic calorimeter (ECL). All these sub-detectors are located 
inside a super-conducting coil that provides a magnetic field of 1.5 T. The return yoke of 
;he coil is instrumented as a K\ and \x detector. A detailed description can be found in 



221. 



This analysis uses 87.7 fb _1 of e + e~ annihilation data taken at the T(4S) resonance 
at ^J~s = 10.58 GeV ( "on-resonance data"), above the production threshold for BB pairs. 
Additional 15.0 fb -1 are taken 60 MeV below the resonance at ^/s = 10.52 GeV ("con- 
tinuum data"), which is also below the production threshold for BB pairs. Hadronic 
events are selected as described in j^J. The selection efficiency of events originating 
from light quarks (d, u and s) passing this hadronic preselection has been estimated to 
be 84.0%, using 9.6 x 10 6 MC events. For c quarks, the efficiency has been determined 
with 6.6 x 10 6 MC events to be 93.0%. The light quark sample contains almost no true 
candidates, reflecting the small rate for gluon splitting into open charm states, i.e. two 
mesons containing c quarks. 

To estimate the efficiency of reconstructing charmed hadrons and to correct for distor- 
tions due to the finite acceptance of the detector, MC samples of e + e~ — > cc events cor- 
responding to a data luminosity of 217 fb -1 (approximately 2 1/2 times the on-resonance 
data), and e + e~ — > qq (q = u, d and s) events corresponding to 18 fb -1 (approximately 
1.2 times the continuum data), have been studied. The MC samples were generated us- 

n 

ing the QQ98 generator |24| employing the Peterson fragmentation function for c quarks 
and were processed through a detailed detector simulation based on GEANT 3.21 j^j]. 
This sample will be referred to as the generic sample. Special samples of several million 
e + e~ — > cc events were generated with the EvtGen generator using the Peterson as 
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well as the Bowler fragmentation functions and were also run through the detector sim- 
ulation. These samples will be referred to as reweighted samples; see Section IVI CI for 
details about the reweighting procedure. For each charmed hadron used in this analysis, 
a sample was generated where that hadron was forced to decay in the same channel as 
later reconstructed. These samples were reconstructed using the same procedures as for 
data. 

A. Particle Identification 

To minimise possible kinematic biases due to tight selection criteria for identified parti- 
cles, only loose cuts on the particle identification of the stable particles have been applied. 
All particles with mean lifetime longer than 100 ps have been called "stable" . Apart from 
reducing a potential kinematic bias, this increased the reconstruction efficiency at the cost 
of introducing more background, especially in the low momentum region. 

In general, the identification for each track was based on one or more likelihood ratios, 
which combined the information from the time-of-flight and Cerenkov counters and the 
energy loss dE/dx in the drift chamber. Pions and kaons were separated by a single 
likelihood ratio £(K)/(£(K) + C(n)). Charged particles were identified as pions if this 
ratio was less than 0.95 and as kaons if this ratio was larger than 0.05. This overlap 
allowed a charged particle to be identified as both a pion and a kaon, potentially resulting 
in identifying a mother (candidate) particle as its own anti-particle (i.e., a D° —* K~n + 
decay could be identified as a D° — ► n~ K + decay), and therefore overestimating the 
number of candidates. As this misidentification was only possible for neutral particles, 
an additional systematic uncertainty has been assigned for the D° and _D*°; see section 
HV|for details. 

For proton identification, similar likelihood ratios were required to fulfil 
£{p)/{£{p) + £(tt)) > 0.6 and £{p)/(£{p) + £{K)) > 0.6. For the vr°, photon candidates 
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with energies above 30 MeV were combined to form a ir° candidate. Under the assump- 
tion that the tt° candidate decayed at the interaction point, it was required to have an 
invariant mass consistent with the 7r° mass. 

The efficiencies e and misidentification probabilities / for tracks from signal candidates 
under these cuts have been measured in data, and are listed in Table HJ in all cases 
except the proton, e > 95% and / < 26%. For kaons and pions the efficiencies and 
misidentification probabilities have been estimated in bins of the particle's momentum 
from D* + and subsequent D°-^K~tt + decays; for protons, A decays have been used. The 
observed momentum spectra in data have been used to derive the listed numbers. 

TABLE I: Typical efficiencies and misidentification probabilities for tracks from signal can- 
didates used in this analysis. The misidentification probabilities listed under ir^ means the 
probability of mis-identifying it as a 



7T ± K± P 

D° (e\f) (96%|26%) (96%|26%) 

D+ (e\f) (96%|12%) (97%|24%) 

Df (e|/) (98%|17%) (9796|21%) 

A+ (e|/) (98%|15%) (97%|21%) (81%|7%) 



In addition to the requirements on the particle identification, all tracks had to be consis- 
tent with coming from the interaction point (IP). For the slow pion from the >D°7r + 
decay, all track quality and particle identification requirements were removed to increase 
the efficiency. 
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B. Reconstruction of charmed Hadrons 



The reconstructed hadron decay chains used in this analysis are the following: 
D° -> K~ir + , D+ -> K-H+1I+, D+ -> 07T+ (0 -> K+K~), A+ -> p + K-ir+, 
D*+ -> D°7T + (D° -> X-7T+), £>*+ -> D+7T (Z>+ -> K-7T+7T+) and 

The inclusion of charge-conjugate modes is implied throughout this paper and natu- 
ral units are used throughout. For all charmed ground state hadrons, candidates whose 
masses were within 50 MeV/c 2 of their respective nominal mass were considered. For the 
intermediate D° and D + coming from the excited states D*° and D* + a mass window of 
15 MeV/c 2 around the nominal masses of the D° and the D + was chosen. Additionally, 
the selection window for the two excited states was tightened to 15 MeV/c 2 around the 
nominal mass difference between the excited meson and the D° or D + . For the inter- 
mediate from the decay, the mass window was chosen to be 7 MeV/c 2 . Multiple 
candidates for each particle and anti-particle were removed by a best candidate selection. 
Most false D*° and D* + candidates were formed from a true D° and a random slow pion. 
Therefore, the slow pions were used to determine the best candidate. For the neutral slow 
pion, the smallest x 2 °f the vertex fit was used. For the charged slow pion, the smallest 
distance to the IP of all hits used in the reconstruction was used. For all other charmed 
mesons, the selection was based on the particle identification of the kaon. In the rare 
case that multiple candidates were formed with the same kaon, the first candidate was 
randomly chosen. 



III. ANALYSIS PROCEDURE 



There are two variables commonly used in the measurements of fragmentation func- 
tions. These are the scaled energy xe = E can didate/E^didate an d the scaled momentum 
x P = |Pcandidatc|/|p^ndWatel, where E^ datc = |p^| datc | = ^s/4 - in 2 , and m H de- 

ll 



TABLE II: The values for the masses or mass differences and branching fractions for all charmed 
hadrons used in this analysis. The masses are used only to shift the mass or the mass difference 
distributions in order to center their peaks near zero, therefore no errors are assigned. The 
branching fractions are taken from 



hadron 



decay mode 



mass [GeV/c 2 



product branching fraction 



D° 

D+ 

Df 

A c + 
D*+ 



D 



*o 



K~TT + 

(f)TT + 
D°TT+ 

D+ir° 

D°7T° 



1.8645 
1.8693 
1.9685 
2.2849 
0.1455 
0.1407 
0.1422 



0.0380 ± 0.0009 
0.092 ± 0.006 
(0.036 ± 0.009) • (0.491 ± 0.006) 
0.050 ±0.013 
(0.677 ± 0.005) • (0.0380 ± 0.0009) 

(0.307 ± 0.005) • (0.092 ± 0.009) 
(0.619 ± 0.029) • (0.0380 ± 0.0009) 



notes the mass of the charmed hadron. For b quarks at higher CMEs, the scaled energy 
xe is often used. In this case, the mass of the B hadron reduces only slightly the allowed 
range at small xg. For charmed hadrons at 10.58 GeV the range of xe is significantly 
reduced; hence xp is prefered and will be used in this analysis. Unless otherwise stated, 
all variables are given in the e + e~ rest frame, taking into account the different beam 
energies for the on-resonance and the continuum samples. 

For various bins in the range from 0.0 to 1.1 in the scaled momentum xp, the signal 
yield has been determined from a fit to the mass or mass difference distributions of all 
candidates within the aforementioned selection windows. The finite momentum resolution 
of the detector can result in events being recorded in the region above the naive limit 
of xp=l, however, in the case of D*, the principal contribution is due to the process 
e+e" -> D*D. See Section EZ3 for details. 
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A bin width of 0.02 in xp has been chosen for all particles as a compromise between 
the statistical precision in each bin and the momentum resolution, which is a factor of 
two smaller. Additionally, to investigate the high xp region around and above the naive 
limit of xp=l, the bin width has been decreased to 0.01; an expanded view of the region 
0.90 < xp < 1.05 with this binning will be discussed in Section IV CI Since this decreased 
bin width is still larger than, but comparable to, the momentum resolution, an unfolding 
using the singular-value-decomposition (SVD) approach (^tJ was tried in addition to the 
normal bin-by-bin correction and is discussed in Section IV CI 

The mass or mass difference distributions were parametrized by a single Gaussian, 
except for the D* + — >D°n + decay channel where a double Gaussian was employed. For the 
mass distributions, the background was parametrized by a quadratic function; for xp > 0.9 
a linear function was found to be sufficient to fit the considerably lower background. For 
the mass difference distributions of the excited D mesons, a phase-space-like function 
/(Am) = a(Am — AMo) b was used with a and b being free parameters and AMq the 
nominal difference between the mass of the excited mother particle and that of the ground 
state charm meson. 

For all charmed hadrons, the mean mass and the width of the signal Gaussian 
was fitted separately for MC, continuum and on-resonance data. For these fits, xp was 
divided into 4 bins from 0.2 < xp < 1.0 with a constant bin size of 0.2. In a second fit, 
two quadratic functions mj(xp) and yUj(xp) were fitted to the results of the first fit in these 
four bins. For the distributions with a bin width of 0.02 and 0.01 in xp, the mean and 
width parameters in the fit were fixed to the values of the quadratic functions rriiixp) and 
/ij(xp) for the appropriate xp value. 

For the D* + ^D°n + decay mode full correlations between the two Gaussians of the 
signal function were taken into account when determining the fit yield. 

When combining the on-resonance data with the continuum data, two corrections have 
been applied to the on-resonance data. After normalising using the integrated luminosities 
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of the respective samples, the naive 1/s dependence on the total hadronic cross section 
has been taken out by multiplying the distributions of the on-resonance sample by the 
square of the ratio of the CME's, namely by (10.58 GeV/10.52 GeV) 2 . Second, from MC 
an additional correction of +0.27% due to different initial state radiation (ISR) at the two 
energy points has been applied to the on-resonance samples. This correction was based 
on a MC study of the total cross sections at these two energy points. 

A. xp-dependent Mass Fits 

Fig. ^ and Fig. |21 show the mass distributions of all charmed hadrons reconstructed in 
this analysis for two representive bins in xp. The xp bins shown are 0.28 < xp < 0.30 in 
Fig.^and 0.68 < xp < 0.70 in Fig- El They represent a low xp bin with higher background 
and a bin close to the maximum of the xp distribution with less background, respectively. 
All mass (mass difference) distributions have been shifted by their nominal mass (mass 
difference) to center the peaks at zero. See table [H] for the masses used. Note that the 
scale on the y-axis does not start at zero in the upper four plots in Fig. ^ and Fig. El 

For 0.28 < xp < 0.30 (shown in Fig. ^), the mass distributions for the D° and D + 
ground states and the mass difference distributions for the excited states show clear peaks 
at the expected value for signal. Compared to higher xp values, the background is higher 
due to a larger amount of combinatorial background, and the signal-to-background ratio 
is lower. At higher x P values, such as those shown in Fig. |21 (0.68 < x P < 0.70), the back- 
ground is considerably reduced, whereas the signal yield is enhanced. This significantly 
increased the signal-to-background ratio. 
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B. Raw Signal Yield 

Fig. El shows the signal yields as a function of xp for all charmed hadrons, not corrected 
for the reconstruction efficiencies and for the branching fractions, denoted with "B" in 
the plots. For all particles, the contribution from B decays is clearly visible in the low xp 
range, which is xp < 0.5 for all charmed mesons containing a light quark as the spectator. 
For Df from B decays, the upper bound is approximately xp ~ 0.4, reflecting the energy 
required to produce an additional strange quark. Contributions from the b — > u transition, 
where the Df is formed at the upper vertex, can populate the region up to x P = 0.5, but 
are strongly suppressed. For the A+, the only baryon reconstructed in this analysis, the 
upper bound is further decreased to approximately x P ~ 0.37, due to the production of 
an additional anti-baryon. 

All distributions peak around xp ~ 0.6 — 0.7 and show similar shapes. 



C. Efficiency Correction 

The efficiencies were determined from MC and are defined as the appropriate raw signal 
yield (determined by the same procedure as for data) divided by the generated MC xp 
distribution. The seven histograms in Fig. 0] show the xp-dependent efficiency of each 
charmed hadron used in this analysis for continuum data and on-resonance data. The D° 
efficiency is close to 50% and almost constant over the entire xp range. The efficiency for 
D* + — >D°7T + approaches the D° efficiency at high values of xp and diminishes at lower 
values of xp, reflecting the reduced efficiency of reconstructing low-momentum pions. 
The two D* decay modes that include a neutral slow pion show a different behaviour: the 
efficiencies stay constant over a wide range of about 0.3 < xp < 1 and below 0.3 < xp the 
efficiency increases for xp— >0 due to the increasing reconstruction efficiency for slow 7r°. 
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The reconstruction efficiencies for the three-particle decay modes do not show a strong 
dependence on xp, slightly varying between 15% and 20% for the and remaining 
constant at about 30% for the A+. The decreasing efficiency for particles at values close 
to the kinematic limit is an artefact of the decreasing statistics in all generic MC samples. 
The reweighted samples, which were generated with a different fragmentation function 
than the generic samples, contain significantly more events in the very high xp region and 
do not show such behaviour. This difference between the two efficiency estimates was 
added to the systematic uncertainty. 

The efficiency is a function of the production angle, which differs for charmed hadrons 
from B decays and from continuum events. For the on-resonance samples, the efficiency 
has been determined by a luminosity-weighted mixture of charmed MC and dedicated 
samples containing decays of charged and neutral B mesons. For the continuum sample, 
only charmed MC was used. 

In data, it was verified that D* + produced in e + e~ annihilation are unpolarised by 
verifying that the distribution of the cosine of the helicity angle is flat. The helicity angle 
is defined as the angle between the slow charged pion in the D* + rest frame and the flight 
direction of the D* + in the center of mass system of the event. Because the efficiency for 
D* + ^D°7i + strongly depends on the momentum distribution of the slow ir + , which in 
turn depends on the helicity angle, polarised D* + can introduce a bias into the efficiency 
correction. 



IV. SYSTEMATIC UNCERTAINTIES 

Various sources of systematic uncertainties have been considered: 
Uncertainties due to tracking were estimated to be 1% per track using a sample of 
partially reconstructed D* + decays. As the uncertainty increased at very low momentum, 
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the estimated momentum-dependent uncertainty of the slow charged pion was folded with 
the observed momentum spectrum. The systematic uncertainty due to the slow neutral 
pion detection efficiency was assessed by examining the differences in the shapes of the 
fragmentation function of the two D* + decay modes, D* + — > D°tt + and D* + — > D + n°. 

Uncertainties due to the modeling of ISR in the MC were determined by restricting the 
longitudinal momentum in the laboratory frame of all candidates to p l ^ b > only. This 
cut preferentially removed events with ISR photons in the negative z direction, potentially 
introducing an artificial asymmetry. The z direction is defined as being anti-parallel to 
the positron beam, which coincides up to corrections due to the crossing angle with the 
boost vector into the e + e~ rest-frame. 

The cut on the likelihood ratios for kaon and proton candidates was tightened to 0.2 and 
0.8, respectively, and the difference was taken into account in the systematic uncertainty. 

Potential differences between the actual signal shape and the fitting function were 
estimated by determining the signal yield with a counting method instead of using the 
fit. Here, the number of entries in the mass (mass difference) distribution was counted in 
a window about one third the size of the total 50 MeV/c 2 (15 MeV/c 2 ) window around 
the peak position, corresponding to roughly three times the resolution. The number 
of background events was subtracted after integrating the background function of the 
standard fit within the same mass window. 

An additional flavour assignment systematic uncertainty was taken into account for 
the neutral states D° and D*°. The loose cuts on the charged pion and kaon particle 
identification allowed a D° to be identified as a D°: the flavour of the D° from D* + 
decays was identified by the charge of the slow pion (except for a small contribution 
from doubly-Cabibbo-suppressed decays). In the MC sample, the likelihood ratio of the 
pion candidate was larger than that of the kaon candidate for 1.3% of all D° candidates; 
the corresponding fraction was determined to be 1.1% for D* + decays. The statistical 
uncertainties on these numbers are less than 0.05%. Accordingly, a difference of 0.2% was 
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assigned as the uncertainty of the flavour assignment due to the overlap of the pion and 
kaon likelihoods of the particle identification. 

The luminosity of the data sample was determined to have an uncertainty of about 
1.4%. A corresponding scale uncertainty of 1.4% was assigned to the normalisation of the 
shape. It has been checked that the normalisation of the fragmentation functions of the 
on-resonance and continuum sample agree with each other, and their difference of 0.94% 
is well within the scale uncertainty. 

Finally, the reconstruction efficiencies of the generic and the reweighted samples dif- 
fered slightly. This small difference was added to the systematic uncertainty. 

All systematic uncertainties were added in quadrature to give the total systematic 
uncertainty. 

V. RESULTS 

In this section, various results for the charmed hadrons are presented. 
A. xp Distributions 

Fig. |S] shows the efficiency-corrected xp distributions for the different particles for e + e~ 
annihilation events, i. e. spectra of hadrons formed in the fragmentation of charm quarks. 
Above xp > 0.5, the differential xp distributions of the on-resonance sample and the 
continuum sample have been combined by a weighted average, where the inverse of the 
squared statistical uncertainty was used as the weight. As the systematic uncertainties 
for both samples are highly correlated, the larger uncertainty of the on-resonance and the 
continuum samples was used for the combined sample. 

As most efficiencies do not depend strongly upon xp, the shapes of the efficiency- 
corrected distributions are similar to those of the uncorrected distributions. All distri- 
butions peak around x P ~ 0.6 — 0.7. To determine the peak position, a direct fit of the 

18 



data to the Peterson fragmentation function was tried. The shape of the data agreed very 
poorly as this model does not include gluon radiation or decays from higher resonances. 
Therefore, a Gaussian function was used to determine the peak position. The fit ranges 
were chosen from xp = 0.4-0.8. The results of the fits are listed in Table ITTT1 together with 
the statistical and systematic uncertainties. The statistical uncertainty was determined 
by the RMS of the distribution divided by y/N. 



TABLE III: The peak positions of all hadrons, fitted with a Gaussian near the peak position. 
The fit range was 0.4 < xp < 0.8. Above xp > 0.5, the continuum sample and the on-resonance 
sample have been combined. The given errors are the statistical and systematic uncertainties, 
respectively. 

„ PEAK 
Xp 

D° -» R-it + 0.587 ±0.001 ±0.002 

D+^K-it + it + 0.600 ±0.001 ±0.001 

D+ -> 0vr+ 0.681 ± 0.002 ± 0.003 

A+^p + K-ir + 0.612 ±0.001 ±0.004 

D* + -> D°7T + 0.631 ±0.001 ±0.002 

D + tt° 0.618 ±0.011 ±0.023 

D*° -> L>°7T 0.631 ± 0.001 ± 0.003 



B. Average Number of Charmed Hadrons per B Decay 

The xp distributions of the on-resonance and continuum samples differ in the con- 
tribution from B decays for xp < 0.5. Fig. shows this difference: the differential xp 
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distribution of the continuum sample was subtracted from that of the on-resonance sam- 
ple. Thus, up to statistical fluctuations it contains only contributions from decays of B 
mesons. 

Table I1VI lists the average number of charmed hadrons per B meson decay together 
with the present world average [28]. In order to determine the average number of charmed 
hadrons produced per B decay, we take the difference between the production rate in the 
on-resonance and the continuum sample and normalise by the B meson production cross 
section, which is estimated to be (1.073 ± 0.019) nb based on the measured luminosity 
and the measured number of BB pairs in this sample. Note that this visible production 
cross-section depends strongly upon the energy spread of the accelerator. The uncer- 
tainties in Table HVI are from the limited statistics (first), the systematics as discussed 
in Section HVl (second), and the luminosity measurement and the uncertainties on the 
branching fractions (third). Note that the luminosity measurement and the determina- 
tion of the number of BB are strongly correlated. Both values agree well within one 
standard deviation with each other, only the average number of produced D* 0, s here is 
lower by about one standard deviation and is closer to that of D* +, s. 

The small bump seen in the xp distributions of the charmed mesons except the at 
xp = 0.35 is due to two body decays of the B mesons such as B — > D°D^*' in case of the 
D°. 



C. High xp Region 

An expanded view of the high xp region is shown in Fig. [7| The downward triangles 
show the efficiency-corrected data; the upward triangles show the corrected and unfolded 
data. 



Unfolding was done using the singular- value-decomposition (SVD) method 27[. From 
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TABLE IV: The average number Nb-+ c of charmed hadrons per B meson decay, corrected for 
acceptance and reconstruction efficiencies. The listed uncertainties are statistical, systematic, 
and the one due to the uncertainties on the branching fractions of the decays involved as well 
as on the luminosity, respectively. 



N B ^c PDG(2004) 



D° — » 


K-TT+ 


0.644 ± 0.003 ± 0.024 ± 0.021 


0.640 ± 0.030 


D+ -» 




0.248 ± 0.004 ± 0.033 ± 0.020 


0.235 ±0.019 




(f>7T + 


0.122 ± 0.015 ± 0.033 ± 0.030 


0.105 ±0.026 


A+- 




0.042 ± 0.011 ± 0.033 ± 0.018 


0.064 ±0.011 


D*° 


D°TT° 


0.217 ± 0.014 ± 0.020 ± 0.018 


0.260 ± 0.027 


D* + 


D°TT+ 


0.218 ±0.007 ±0.020 ±0.015 


0.225 ± 0.025 




D + TT° 


0.202 ± 0.014 ± 0.022 ± 0.018 


0.225 ± 0.025 


average 


D*+ 


0.215 ± 0.006 ± 0.022 ± 0.016 


0.225 ± 0.025 



MC, we determined the response matrix of the detector for producing for a certain true 

input value of xp true a measured value of x P measure d- This matrix was decomposed using 

the SVD into two orthogonal and one diagonal matrices which can easily be inverted. 

I I 

Inverting the diagonal matrix was limited by a criteria denned in [27] to contain only 
elements, which are of statistical significance. 

The hatched histogram show the only process e + e~ — > D* + D~, the open histogram 
shows the sum of the previous process and e + e~ — > D* + D^~ . 

The xp distributions for the ground states D°, D + , and extend up to the naive 
kinematic endpoint xp = 1 and no significant number of events are present for x P > 1. 

All three xp distributions for the excited D mesons, however, show an enhance- 
ment at xp > 1. These events above xp = 1 correspond to events of the processes 
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e + e — > D* + D^ or e + e — > D*°D^*>° and are in good agreement with the measured 
cross sections Q of 0.55 ± 0.03 ± 0.05 pb for e + e~ -»■ and 0.62 ± 0.03 ± 0.06 pb 

for e + e~ — > D* + D*~ . Note that these events populate xp > 1 only because of the use 
of the simplified upper limit |p^ndidatel' f° r producing two D* mesons. A background 
fluctuation producing an artificial peak is unlikely for two reasons. First, at high xp, 
the background is negligible, and second, the unfolding procedure tends to identify sig- 
nals at the edge of a distribution as statistical fluctuations rather than real signals, thus 
decreasing the significance of the signals. 



D. Total Production Cross-Section 

The total production cross-section is given by the integral of the xp distribution. This 
integral was determined for the continuum sample using the current value of the world 
average product branching fraction of each particle, see Table HT1 and [28]. The results are 
listed in Table W\ where the third error component reflects the uncertainty on the product 
branching fraction. 

The results by CLEO Q g ,ven in the last coln m n used the, owe branching .actions, 
which differ slightly from the world averages used here. The results, however, agree 



30j is given in the same column. 



well with each other. Another measurement by BaBar 
The total production cross-section for the D*° differs only slightly from that of the D* + . 
This can be understood as resulting from different feed-down contributions from higher 
resonances. 
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TABLE V: The total production cross-sections e + e~ — ► DX (or A+X), which have been cor- 
rected using the current world average of the respective product branching fractions. The listed 
uncertainties are statistical, systematic and the uncertainty due to the uncertainty on the branch- 
ing ratios. Other measurements of production cross-sections are listed. The third column shows 
measurements by CLEO 

Q, BaBar Q (marked «) and an older CLEO measurement Q 

(marked ( 2 **). 





X 


CPROD [pb] 


a PROD(CLEO'04/ BaBar) [pb] 


D° — > 


K--K+ 


1449 ± 2 ± 64 ± 38 


1521 ± 16 ± 62 ± 36 


D+ -» 


K~TT + TT + 


654 ± 1 ± 36 ± 46 


640 ± 14 ± 35 ± 43 




(f>7T + 


231 ± 2 ± 92 ± 77 


210±6±9±52( 1 ) 






189 ± 1 ± 66 ± 66 


270 ± 90 ± 70( 2 ) 


D *o _^ 


D°TT° 


510 ± 3 ± 84 ± 39 


559 ± 24 ± 35 ± 39 


D* + 


D°TT+ 


598 ± 2 ± 77 ± 20 


583 ± 8 ± 33 ± 14 


L>*+ 


D + TT° 


590 ± 5 ± 78 ± 53 




average 


D*+ 


597 ± 2 ± 78 ± 25 





E. Mean Values for xp and Moments 

In addition to the peak position for the seven xp distributions, the mean and higher 
moments of these distributions were determined from distributions in (xp) n with a bin 
width of 0.02 in xp and a bin-by-bin efficiency correction was applied. The n th moment 
was determined by the mean of the efficiency corrected distributions in (xp) ra , and its 
statistical uncertainty was determined by a/y/1%, where is the number of entries in 
the uncorrected (xp) n distribution. Tables fVTl and fVTTl show the moments for the different 
decay modes. 
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TABLE VI: The first three moments of the xp distribution for the seven particles/decay modes 
used in this analysis. The listed uncertainties are statistical and systematic uncertainties. 



D° -> K- vr+ 

D+ -> K- 7T+ 7T+ 

D+ -> <f> 7T+ 

A+ -> p + if" 7T+ 



(x P ) x 1000 

570.33 ±0.18 ±2.23 
578.03 ± 0.18 ± 1.47 

635.34 ±0.47 ±4.22 
582.45 ± 0.39 ± 2.53 
612.17 ± 0.36 ± 1.43 
586.06 ±0.37 ±16.10 
607.63 ± 0.42 ± 6.07 



(x P 2 ) x 1000 

353.98 ±0.29 ±2.50 
363.42 ±0.29 ±1.58 
442.52 ±0.81 ±8.64 
364.94 ± 0.63 ± 3.40 
407.96 ±0.61 ± 2.01 

380.99 ± 0.64 ± 17.89 
401.98 ±0.69 ±5.60 



(x P 3 ) x 1000 
233.85 ± 0.28 ± 2.54 

243.58 ± 0.27 ± 1.58 
323.52 ±0.83 ±11.14 

239.59 ± 0.59 ± 2.24 
286.97 ±0.60 ±3.38 
266.49 ± 0.62 ± 17.05 
282.65 ± 0.68 ± 5.90 



F. Production Angle 

Taking the interference between the exchange of virtual photons and Z bosons into 
account, the differential cross section for e + e~ — > cc is modified from the 1 + cos 2 6 form 
of the Born amplitude for pure photon exchange: 

<]<1 |(1 + cos 2 6)a T + - sin 2 0a L + - cos 8a A (1) 



dcos9 8 V ' 4 4 

Here, 6 describes the angle between the incoming electron beam and the outgoing hadron 
containing the charmed quark, as measured in the CM frame. The term a? describes 
the contribution of pair production of spin-1/2 particles from transverse polarised vector 
bosons, the term ctl the contribution from longitudinal polarised vector bosons and the 
term a a denotes the parity violating asymmetry due to the interference between Z bosons 
and virtual photons. 

The /C/C MC generator 0| was used to predict the production angle distributions for 
the different charmed hadrons. This MC generator includes interference between initial 
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TABLE VII: The fourth through the sixth moments of the xp distribution for the seven parti- 
cles/decay modes used in this analysis. The listed uncertainties are statistical and systematic 
uncertainties. 



(x P 4 ) x 1000 



(x P 5 ) x 1000 



(x P 6 ) x 1000 



D° -» K--K+ 

D+ 

A+ p+A'-vr 4 



161.83 ±0.25 ±2.19 

171.54 ± 0.24 ± 1.28 
244.69 ± 0.81 ± 12.06 
163.04 ±0.52 ±5.11 

211.55 ±0.57 ±5.36 
196.24 ± 0.58 ± 16.67 
215.63 ± 0.68 ± 12.47 



116.97 ±0.22 ±1.90 
125.62 ±0.22 ±1.16 
188.72 ±0.76 ±10.64 
115.07 ± 0.46 ± 1.94 
162.26 ±0.53 ±7.05 
150.28 ± 0.53 ± 15.30 
160.05 ± 0.59 ± 7.98 



88.08 ±0.19 ±4.77 
95.52 ±0.19 ±1.11 
150.59 ±0.72 ±9.75 
85.06 ±0.41 ±2.31 
128.24 ±0.49 ±8.17 
118.85 ±0.49 ±13.89 
126.87 ±0.54 ±8.72 



and final state radiation (ISR and FSR) as well as electro-weak corrections. 10 8 cc events 
were generated with /C/C and hadronised with the PYTHIA generator. 

For the generated events, x P was divided into 20 bins of equal width and a three- 
parameter fit to the production angle was performed: 



f(0,x P ) 



a (x P j 



3 3 3 

-(1 + cos 2 ^) + -r s (xp) sin 2 # + -r c (xp) cos9 

8 4 4 



(2) 



where ao is the normalisation and r$ and rc are the relative contributions for the sin 2 9 
and the cos 9 terms, respectively. The results of these fits are shown in Fig. |S] as the solid 
line, together with the measured data points. 



For data, the signal yield was determined in bins of xp and cos 6, where 9 is the 
production angle of the charmed hadron. It should be noted that here, the efficiency 
correction depends on xp and cos 9. 
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xp was divided into only 4 bins: xp < 0.3, 0.3 < xp < 0.5, 0.5 < xp < 0.7 and 0.7 < xp. 
The boundaries were chosen in order to roughly equalise the number of candidates per 
bin. Two bins each were chosen below and above xp = 0.5, which is the upper kinematic 
limit for hadrons from B decays. 

cos9 was divided into 20 bins. In each bin of xp and cos#, the efficiency corrected 
signal yield in the mass or mass difference distributions was fitted separately for the on- 
resonance and continuum samples. The same signal and background functions were used 
as in the fits which depended only on xp. 

In a second step, a three-parameter fit (similar to Eq. EJ) to the signal yields was 
performed in bins of xp. The fit range was restricted to —0.8 < cos# < 0.8. As a 
systematic check, the range was tightened to —0.7 < cos9 < 0.7. 

No significant deviation of the r$ and the rc parameters from their expectation was 
found for the continuum sample. The expectation from the /C/C generator was of the 
same order as the statistical uncertainties except for the rs term in the lowest xp bin, 
where gluon radiation introduces a longitudinal momentum component and therefore 
smears out the initial distribution of the production angle. This smearing introduces a 
significant sin 2 6 term, to which the measured values in this regime agree well. As the 
number of entries in the low xp bins also diminish, the statistical uncertainties increase 
to roughly the same size as the expected effect. The fitted values for rs and rc in the 
lowest xp bin in the continuum samples suffer from very low statistics (see Fig. 0}. For 
the Df and the A+, the lowest bins in the continuum samples have been neglected as the 
numbers of entries in these bins were too low to perform an angular analysis. 

For the on- resonance sample with higher statistics, the rs terms significantly deviate 
from zero for xp < 0.5, as B decays with different production angle distributions also 
contribute. The rc term is again consistent with both zero and the expectation from 
the MC generator. Tables IVIIII and llXI list the measured values for the rs and rc values 
for the continuum and the on-resonance data. The systematic uncertainties here include 
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the uncertainties discussed in the standard analysis as well as the uncertainty due to the 
restricted fit range. 

VI. INTERPRETATION 

A. Contributions from Higher Resonances 

Contributions from excited states have been considered only in the xp distributions 
of the D° and the D + and for these, only contributions from D*° and D* + were consid- 
ered. For higher resonances, both production cross-section and branching fractions of e.g. 
D** have large uncertainties and have been neglected. In order to reduce the statistical 
uncertainty, a MC-based correction was applied: Three large samples of several million 
MC events were generated without a detailed detector simulation. These samples were 
required to contain the decay modes D* + ^D°tt + , D* + ^D + ti and D*°— >D°7r°, respec- 
tively. For these events, the xp of the parent D* + /D*° vs. the xp of the daughter D°/D + 
were stored in a two dimensional matrix. The measured and efficiency corrected D*° 
and D* + xp distributions were multiplied with this matrix in order to estimate the xp 
distribution for all D 0, s and D +, s coming from D* + /D*° decays. 

The two plots at the top of Fig. |U] show the contributions of D* + and D*° decays 
to the D° fragmentation function (left), and the contribution of D* + decays to the D + 
fragmentation function (right). These plots are not corrected for the branching fraction 
of the D decay. The bottom plot in Fig. |U] shows primary D° and D + fragmentation 
spectra: the total D° (D + ) spectrum minus the contribution from D*° and D* + (D* + only) 
decays. The difference of 13% between the sum of primary produced D° and D + should be 
compared to the 6.5% relative uncertainty in the D + ^K~h + 'k + branching fraction. Also, 
as only the contribution from D* decays has been considered, the remaining difference 
may be due to the contribution of other resonances. In the generic MC sample, where 
primary _D°'s and D + 's were produced in equal amounts, there was an excess of 6% in the 
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production of D° (compared to D + ) mesons in the decay of resonances other than D* + . 



B. Ratios 

Comparisons of production rates for various particles are useful for understanding the 
dynamics of fragmentation, as systematic errors cancel in the ratio. In this section we 
present ratios of both integrated cross-sections and cross-sections as a function of xp, to 
characterise general properties of fragmentation and to test the agreement between MC 
simulation and data. 

Table |Xj presents three ratios of total production cross sections. Since the production 
of D* + and D*° is included in the total D + and D° production rate (all D* + decay to 
either D + or D°, and all D*° to D°), the ratio of D* to D production measures V/ (V + P), 
the probability of producing a vector charmed meson. (Here we write V for vector and P 
for pseudo-scalar meson production rates.) A correction is necessary to account for higher 
resonances decaying directly to D +, °. For example, based on the measured production 

n n 

rates of the _Di(2420) and 1^2(2460) [32] mesons, known branching fractions [28], and 
isospin relations, we find a correction of —(3.7 ± 3.3)% to the first ratio. In principle, 
further corrections due to decays of broad D** states and charmed-strange mesons are also 
required. However, no corrections have been applied to the values presented in Table 1X1 

Similarly, the second ratio measures the production rate of charmed-strange mesons as 
a fraction of all charmed mesons, up to corrections for D s i(2536) and -D* 2 (2573) decays. 
The third ratio measures the production rate of charmed baryons relative to that of 
charmed mesons, excluding the charmed-strange baryon states. For comparison, see j^j. 

Ratios of production rates as a function of xp allow momentum-dependent effects in 
fragmentation to be studied, although contributions from decays of higher states also 
appear. Fig. 1101 shows the following five ratios as a function of xp, for both on- resonance 
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and continuum data: 

(a) x P (D* + )/xp(Dp rim ), sensitive to the production rate of vector relative to pseudo- 
scalar mesons; 

(b) x P (D° im )/xp(D^ rim ), sensitive to charged relative to neutral pseudo-scalar produc- 
tion; 

(c) xp(D+)/xp(Dp rim ), sensitive to the production of strange quarks; 

(d) xp(A+)/xp(Z)i.j m ), sensitive to the production of baryons relative to mesons; 

(e) xp(D*°) /xp(D* + ) , the relative production rate of the vector mesons. 

The suffix "prim" denotes xp distributions corrected for the contributions from D* decays; 
D + production has been chosen as the denominator in (a)-(d), as this correction is smaller 
than that for D°. No other corrections for excited states have been applied. 

The production ratios in Fig. ITUT a) and (b) are similar for on- resonance and continuum 
data. In Fig. HUf c). the contribution of B meson decays to production can be clearly 
seen in the low-xp region. In Fig. ITOT d). baryon production in B decays is seen to be 
suppressed below xp w 0.4. As xp approaches unity, the A+/D + production ratio goes to 
zero, consistent with the conservation of baryon number. 

Four similar ratios are shown in Fig. ITTT aWd) for both continuum data (full squares) 
and MC simulations, to test the performance of the MC for various fragmentation function 
parameters. In these plots, the total D + production rate, without D* subtraction, is used 
in the denominator of the ratios. The open histograms show the generic MC sample, which 
agrees with the data only for the highest values of xp of the distributions in Fig. UTT a) and 
(b), but fails to describe the data distributions at lower values. The open squares show 
a second MC sample generated with the Bowler fragmentation function, which shows a 
similar behaviour. 
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Noting that the parameter PARJ(13) in PYTHIA gives the probability for a charmed 
hadron produced in fragmentation to have spin one, 50 MC samples of 10 7 events each 
were generated, with PARJ(13) values ranging from 0.25 to the default value 0.75 given 
by spin counting. These samples were generated in addition to the "reweighted samples" 
used for more refined MC comparisons as described in the next section. A reduced chi- 
squared x 2 was calculated for these samples and the measured and corrected ratios, and 
a fourth-order polynomial in PARJ(13) was fitted to the results. The minimum y 2 
was found to occur at PARJ(13) = 0.592 ± 0.021 for the x P (D* + ) /x P (D + ) ratio, and 
0.592±0.046 for the xp(D°) /xp(D + ) ratio, where the uncertainties denote the la range in 
the fitted polynomial. We note that models of hadron production more sophisticated than 
spin counting predict values for PARJ(13) below 0.75; see 34] and references therein. 

In Fig. II If a Wd) a third MC sample generated with the Bowler fragmentation function, 
and PARJ(13) = 0.59, is shown with closed triangles. This sample and the data agree 
well within the error bars over almost the entire range in Fig. lTlT a) and (b). All three MC 
samples fail to describe the ratios in Fig. ITIT c) and (d): both the endpoints and the shape 
disagree. The difference between the MC samples is small compared to the discrepancy 
with data for Df / D + production in Fig. [TTT c); while Bowler fragmentation (open squares 
and triangles) gives an improved description of A+/D + production in Fig. ITTY d). the 
agreement is still poor. There are no obvious parameters in the MC which can affect 
these two ratios in such a way as to improve the agreement between data and MC. 



C. Comparison of xp distributions with predictions from MC generators 

The models used by MC generators are based on simplified assumptions and require 
input from experiment: this is reflected in the models' input parameters. The commonly 
used JETSET/PYTHIA generators are based on the Lund or string model, in which a 
coloured string is expanded between two emerging partons. The energy stored in the 
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string increases with increasing distance, and eventually allows the production of a new 
quark anti-quark pair. The quark (anti-quark) then produces a meson together with the 
initial anti-quark (quark). The energy distribution of the new quark or anti-quark is 
described by a fragmentation function. Various fragmentation function models have been 
published; see Table IXT1 for a summary. These models depend upon up to two independent 



variables, these are the transverse mass m± = y m 2 + p 2 ± of the newly created hadron, 
and z, the fraction of the longitudinal energy E + p« which the meson inherits from the 
initial quark. 

Not all models listed in Table IXll are implemented in the JETSET/PYTHIA gener- 
ator. In order to be able to compare all models to data, a reweighting technique has 
been applied. Here, several million events referred to as "reweighted samples" have been 
generated, allowing a more elaborate comparison than described in Section IVI Bl For 
these events, z and p± were stored together with the event. This allowed each event to 
be reweighted in order to mimic any other fragmentation function. Scans through the 
parameter space of the five listed fragmentation functions have been performed on these 
special samples. This analysis was performed on five different hadrons; D*°— >D°7r° and 
D* + ^D + 7T° have been omitted because of the large systematic uncertainty due to the 
detection efficiency of the slow neutral pion. 



For data and MC, the xp distributions were compared using uncorrected data and the 
reweighted special MC samples after full detector simulation. A x 2 was calculated based 
upon the distribution of the reweighted special sample and the measured data distribution. 
Only statistical uncertainties in each xp bin were taken into account and only bins which 
contained entries in data or MC were included. The number of bins minus the number of 
parameters of the fragmentation function was used as the number of degrees of freedom 
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(d.o.f.). The weights in the reweighting procedure were constructed in such a way that 
the number of events before and after reweighting stayed constant. This way the total 
value of the % 2 becomes dependent on the size of the data and MC samples; the relative x 2 
values, however, allow a direct comparison between the different fragmentation functions. 

Table IXTTl shows the Xmin/ d-o.f. for all five particles and five fragmentation functions. 
For all five particles a similar trend is visible. The Bowler model in general agrees best 
with the data. The Lund models shows a similar performance in describing the spectra, 
its Xmin/d-o.f. being by factors of 2-3 better than the next best model. For D + and 
D* + the Xmin/d-o.f. is slightly worse than for the Bowler model. In the minimum of the 
X 2 / d.o.f. distributions the a parameter deviates strongly from the default for most of the 
particles. As the Lund model is employed for fragmentation of all flavour species, such a 
large change in the parameter would also change the particle spectrum of light mesons. 
Therefore, further tuning of the second parameter to the Lund fragmentation function 
has been omitted. 

The models by Collins and Spiller and by Kartvelishvili show a similar Xmin f° r au 
particles, about factors of two to three worse than that of the best models. The last 
model, that of Peterson, shows the worst agreement with a reduced Xmin °f 15 and well 
above, ruling out this model for describing data at this CME. 

The input parameters for the fragmentation functions at the minimum of the x 2 /d.o.f . 
distributions are listed in Table 1X1111 

In summary, the Bowler model shows the best agreement between data and MC, how- 
ever, large differences are still present. These differences might be resolved by adjusting 
other parameters of the generators as well, but such a task is out of scope for this analy- 
sis. The Lund model shows the second best agreement. The models by Kartvelishvili and 
by Collins and Spiller show larger deviations and the commonly used model by Peterson 
shows the worst agreement between data and MC. 
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VII. SUMMARY 



A new determination of the charm fragmentation function at a CME close to the T(4S) 
resonance has been presented. The measured xp spectra have been compared to those 
of five different parametrisations in MC via a reweighting procedure, and the best input 
parameters have been found. The best agreement between data and MC has been found 
for the Bowler model and the Lund model. Additionally, the peak positions and the first 
six moments of the xp distributions have been measured. These measurements will allow 
detailed comparisons between experiment and theory. The total production cross-section, 
as well as xp dependent ratios of the fragmentation functions, place stringent tests on 
existing MC generators, which so far completely fail to describe the xp dependent ratios 
of xp(Df) /xp(D + ) and xp(A + )/xp(.D + ). For the first time, the production rates of D + 
and D° excluding the decay of D* mesons have been measured. They were found to agree 
reasonably well with each other. 

The efficiency corrected data points will be made available via download in the Durham 
HEP REACTION DATA DataBase 0. It is presented in a different way as shown 
this article. Separate sets of the continuum and the on-resonance samples are given as 
sBdcr/dxp, i.e. scaled by the nominal center-of-mass energies of 10.52 GeV and 10.58 
GeV, respectively, and not corrected for the branching ratios. The on-resonance data 
includes the additional correction for ISR of +0.27%, see Section ITTTl for details. 
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FIG. 1: Mass and mass difference distributions for all charmed hadrons reconstructed in this 
analysis, for 0.28 < xp < 0.30 for the continuum sample. The histograms show the data, the 
dotted line describes only the background, the full line includes the signal. The top row shows 
the D° (left) and the D + (right), the middle shows the Df (left) and the A+ (right), and the 
bottom row shows D* + — ► D°tt + (left), the alternative decay mode D* + — > D + 7r° (middle), and 
the D*° (right). 37 
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FIG. 2: Mass and mass difference distributions for all charmed hadrons reconstructed in this 
analysis, for 0.68 < xp < 0.70. The order of the plots is the same as in Fig.^ As in the previous 
figure, the histograms show the data, the dotted lines describe only the background, the full 
lines include the signal. 
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FIG. 3: The signal yield not corrected for efficiencies for the charmed hadrons. The order of the 
particles is the same as in Fig. ^ The contribution from B decays in the on-resonance samples 
(down-left hatching) is clearly visible in the region xp < 0.5. The error bars show the statistical 
uncertainties only. 
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FIG. 4: The efficiencies for the charmed hadrons used in this analysis. The order of the particles 
is the same as in Fig.^ The different production angle distributions for the on-resonance (down- 
left hatching) and the continuum sample (down-right hatching) result in different efficiencies for 
these samples. The error bars show the statistical uncertainties only. 
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FIG. 5: Efficiency corrected momentum distributions for the charmed hadrons produced in 
e + e~-annihilation events, i.e. from fragmentation of charm quarks. The order of the particles is 
the same as in Fig. ^ For xp > 0.5, the on-resonance and continuum data have been combined 
by a weighted average. The inner error bars show the statistical, the outer error bars the total 
uncertainties. 
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FIG. 7: An expanded view of the high xp region. The downward (upward) triangles show 
the efficiency-corrected (unfolded) histograms. Events above the naive limit of xp= 1 can be 
produced via e + e~ — * D* + D^; this is shown in the D* + xp distributions as a hatched (open) 
histogram for e + e _ — * D* + D~ (sum of e + e~ — > D* + D~ and e + e _ — * D* + D*~). The order of 
the particles is the same as in Fig. ^ The inner (outer) error bars show the statistical (total) 
uncertainties. 



43 



0.8; 
0.6 E 
0.4 E 
0.2 E 

°E 
-0.2 E 
' 0.2! 
0.151 
0.1 i 
0.05! 

-0! 
-0.05! 

-0.1! 
-0.151 





D° 








*i 


ft tt 




* — ■ — 







0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x p 





it 



a; 



' 0.2- 
0.1 5 h 
0.1 r 
0.05 r 

-0!f-* 
-0.05 h 
-0.1 r 
-0.151 

0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 



1 

0.8 
0.6 
0.4 
0.2 


-0.2 
' 0.2 
0.15 
0.1 
0.05 

-0.05 
-0.1 
-0.15 




r 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 



1 

0.8 
0.6 
0.4 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x p 



0E 



-0.2 - 
' 0.2? 
0.15! 

0.1! 
0.05! 

-0! 
-0.05! 

-0.1! 
-0.151 





i 1 

— * .i 










' 1' P , 



0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 

x p 



FIG. 8: The production angle coefficients r$ (upper distributions) and tq (lower distributions) 
for the four ground state charmed hadrons: D° and D + at the top, Df and A+ in the middle. 
The sum of all ground states after efficiency correction and an excited D state (D* + — > D°tt + ) 
are shown at the bottom. The upwards (downwards) triangles show on-resonance (continuum) 
data. The connected points show the results from the /C/C generator. For display purposes, the 
points are slightly separated in xp. 
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TABLE VIII: The coefficients in front of the sine-squared term (rs) and the cosine term (rc) 
for different xp bins in the continuum sample. 



particle 



range m xp 



r s 



D° 



0.0 - 0.3 

0.3 - 0.5 

0.5 - 0.7 

0.7 - 1.0 



0.393 ±0.178 ±0.166 
0.035 ±0.016 ±0.007 
0.025 ±0.008 ±0.014 
0.009 ±0.010 ±0.004 



0.048 ± 0.079 ± 0.008 
0.020 ±0.011 ±0.003 
0.010 ±0.006 ±0.001 
0.013 ±0.007 ±0.009 



D+ 



0.0 - 0.3 

0.3 - 0.5 

0.5 - 0.7 

0.7- 1.0 



0.507 ±0.443 ±0.349 
0.039 ±0.034 ±0.008 
-0.040 ±0.011 ±0.006 
-0.003 ±0.012 ±0.019 



-0.076 ± 0.222 ± 0.005 
0.008 ± 0.024 ± 0.002 
-0.028 ± 0.008 ± 0.004 
-0.042 ± 0.008 ± 0.010 



Df 



0.3 - 0.5 
0.5 - 0.7 
0.7 - 1.0 



0.063 ±0.120 ±0.013 
0.076 ± 0.041 ± 0.004 
0.096 ±0.038 ±0.027 



0.073 ± 0.077 ± 0.029 
0.069 ± 0.026 ± 0.007 
0.107 ±0.024 ±0.017 



A+ 



0.3 - 0.5 
0.5 - 0.7 
0.7 - 1.0 



0.072 ±0.112 ±0.027 
-0.023 ± 0.029 ± 0.042 
0.016 ± 0.048 ± 0.037 



0.104 ± 0.077 ± 0.066 
0.006 ± 0.021 ± 0.004 
-0.059 ±0.031 ±0.010 



Sum 
of 
ground 
states 



0.0 - 0.3 

0.3 - 0.5 

0.5 - 0.7 

0.7 - 1.0 



0.654 ±0.312 ±0.142 
0.042 ±0.019 ±0.007 
-0.007 ±0.007 ±0.011 
0.010 ± 0.008 ± 0.007 



-0.019 ±0.124 ±0.039 
0.025 ±0.013 ±0.008 
-0.006 ± 0.005 ± 0.002 
-0.007 ±0.005 ±0.008 



D* 



0.0 - 0.3 

0.3 - 0.5 

0.5 - 0.7 

0.7 - 1.0 



0.221 ±0.157 ±0.779 
0.051 ±0.031 ±0.035 
0.011 ±0.014 ±0.013 
0.011 ±0.015 ±0.012 



0.069 ± 0.070 ± 0.333 
-0.007 ±0.020 ±0.017 
-0.034 ± 0.010 ± 0.005 
-0.056 ± 0.010 ± 0.004 
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TABLE IX: The same as in Table IVTTT1 but now for the on-resonance sample. 



particle 



range m xp 



r,s 



rc 



D ( 



0.0 - 0.3 
0.3 - 0.5 
0.5 - 0.7 
0.7 - 1.0 



0.318 ±0.023 ±0.022 
0.191 ±0.007 ±0.016 
0.005 ± 0.004 ± 0.004 
0.010 ±0.005 ±0.010 



-0.013 ±0.011 ±0.008 
0.001 ± 0.004 ± 0.004 
0.006 ± 0.003 ± 0.002 
0.023 ± 0.004 ± 0.006 



0.0 - 0.3 

0.3 - 0.5 

0.5 - 0.7 

0.7 - 1.0 



0.377 ±0.063 ±0.026 
0.202 ±0.015 ±0.008 
-0.011 ±0.006 ±0.005 
-0.003 ± 0.006 ± 0.007 



0.093 ±0.027 ±0.016 
0.012 ±0.008 ±0.001 
-0.002 ± 0.004 ± 0.001 
-0.014 ± 0.004 ± 0.001 



0.0 - 0.3 
0.3 - 0.5 
0.5 - 0.7 
0.7 - 1.0 



0.428 ±0.130 ±0.032 
0.284 ±0.038 ±0.011 
0.025 ±0.017 ±0.008 
0.026 ±0.017 ±0.008 



-0.079 ± 0.057 ± 0.030 
-0.024 ±0.019 ±0.005 
0.005 ±0.011 ±0.001 
-0.013 ±0.011 ±0.008 



Aj 



0.0 - 0.3 
0.3 - 0.5 
0.5 - 0.7 
0.7 - 1.0 



0.336 ±0.172 ±0.044 
0.048 ± 0.049 ± 0.027 
-0.030 ±0.014 ±0.015 
-0.037 ±0.020 ±0.002 



-0.058 ± 0.084 ± 0.009 
-0.021 ± 0.033 ± 0.020 
-0.023 ±0.010 ±0.003 
-0.024 ± 0.014 ± 0.004 



Sum 
of 
ground 
states 



0.0 - 0.3 

0.3 - 0.5 

0.5 - 0.7 

0.7 - 1.0 



0.349 ± 0.030 ± 0.001 
0.197 ±0.008 ±0.013 
-0.004 ±0.004 ±0.001 
0.002 ± 0.004 ± 0.006 



0.014 ±0.014 ±0.012 
0.004 ± 0.005 ± 0.001 
-0.000 ± 0.002 ± 0.001 
0.001 ± 0.003 ± 0.003 



D 



*+ 



0.0 - 0.3 
0.3 - 0.5 
0.5 - 0.7 
0.7 - 1.0 



0.487 ± 0.067 ± 0.049 
0.401 ± 0.021 ± 0.004 
0.037 ± 0.008 ± 0.002 
0.014 ± 0.008 ± 0.001 



0.014 ±0.030 ±0.012 
-0.004 ±0.010 ±0.008 
-0.030 ± 0.005 ± 0.001 
-0.037 ± 0.005 ± 0.004 
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FIG. 9: Upper plots: The contributions from D* + (up-right hatching) and D*° (up-left hatch- 
ing) decays to the D° (left) and D + (right) xp distributions. These plots are not corrected for 
the branching fraction of the D decay. Lower plot: The xp distributions for primarily produced 
D (error bars) and D + (up-right hatching) mesons. Note that only contributions from D* + 
and D*° decays have been considered; the exclusion of higher resonance contributions can par- 
tially account for the differences in the xp distributions of the D° and D + . Only statistical 
uncertainties are shown. 
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TABLE X: Three ratios of total cross sections, each of the form a(e + e~ — > AX) / a{e + e~ — > 

for the continuum sample. The denominators of all ratios contain the contribution from D* —* D 

and other decays. For comparison, see [3^ 

A B ratio 

D*° + D*+ D+ + L>° 0.527 ± 0.013 ± 0.024 

D+ D+ + D+ + L>° 0.099 ± 0.003 ± 0.002 

A+ Df + D + + D° 0.081 ± 0.002 ± 0.003 



TABLE XI: The functional form of the fragmentation functions used in this analysis. The nor- 
malisation N is different for all functions. The models by Collins and Spiller and by Kartvelishvili 
are not included in the JETSET/PYTHIA generator. 



fragm. function 



functional form 



comment 



Bowler 
Lund 

Kartvelishvili 

Collins-Spiller 

Peterson 



^^(l-^) a exp(-^ 



A±(l-z) a exp(-^i 



a, b identical for all quarks 
a, b identical for all quarks 



Az Qc (l - z) 



1-2 



(l+- 2 )(l-i-A 



Z \ Z 1—2 



widely used 
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FIG. 10: The ratios x P (L>*+)/x P (D+ rim ) and x P (Z^ m )/x P (£>+ im ) at the top, 
x P (L>+)/xp(D+. m ) and x P (A+)/x P (L>+ -J in the middle and x P (L>* )/x P (L>*+) at the bot- 
tom. The open upward triangles represent on-resonance data, and the full squares with error 
bars represent continuum data. The inner (outer) error bars show the statistical (total) uncer- 
tainties. 
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FIG. 11: The ratios x P (D* + )/x P (D + ) in a), x P (D°)/x P (D + ) in b), x P (D+)/x P (D+) in c) and 
xp(A+)/xp(L> + ) in d). The hatched histograms represent the ratio for continuum data, the 
open histogram shows the distribution for the corresponding ratios from generic MC, the open 
squares show the predictions using the Bowler fragmentation function with default parameters. 
The full triangles show the predictions using the Bowler fragmentation function with a tuned 
value for the probability of producing a charmed meson with spin=l: PARJ(13) = 0.59 . 
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TABLE XII: The minimum of the chi-squared distribution, Xmim f° r ^C sam pl es reweighted to 
represent the fragmentation functions shown, varying their respective parameters. The number 
of degrees of freedom (d.o.f.) is also shown for each case. 





D° 


D+ 






D*+ 




xLJd-o.f. 


X 2 min/ d -0-f- 


Xmin/d-O.f. 


Xmvn/d.O.f. 


X 2 min/d-0.f. 


Bowler 


1327.0 / 59 


188.4 / 60 


730.7 / 55 


269.1 / 60 


541.8 / 55 


Lund 


1500.5 / 59 


527.1 / 60 


513.2 / 55 


266.6 / 60 


965.6 / 55 


Collins and Spiller 


3032.1 / 58 


948.0 / 60 


1412.5 / 55 


2836.7 / 59 


1540.7 / 54 


Kartvelishvili 


3210.4 / 59 


861.4 / 60 


735.3 / 55 


390.7 / 60 


1271.1 / 54 


Peterson 


5070.2 / 59 


2229.6 / 60 


829.6 / 55 


1345.0 / 59 


3003.0 / 54 



TABLE XIII: The parameters of the fragmentation functions at the minimum of the \ 2 /d.o.f. 
distributions . 



D° D+ D+ A+ D*+ 



parameter pars at min. pars at min. pars at min. pars at min. pars at min. 



Bowler 


a\b 


0.12 | 0.74 


0.12 | 0.58 


0.12 | 0.68 


0.34 | 0.74 


0.22 | 0.56 


Lund 


a 


0.26 


0.45 


0.2 


0.55 


0.58 


Collins and Spiller 


e' c 


0.04 


0.055 


0.04 


0.04 


0.075 


Kartvelishvili 


a c 


4.6 


4 


5.6 


3.6 


5.6 


Peterson 




0.028 


0.039 


0.008 


0.011 


0.054 
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